A Supersymmetric Solution to the Solar and Atmospheric Neutrino Problems 



O 

o 
o 

(N 



J. C. Romao \ M.A. Di'az 2, M. Hirsch 3, W. Porod 3, and J. W. F Valle ^ 
^ Departamento de Fisica, Instituto Superior Tecnico 
A. Rovisco Pais, 1049-001 Lisboa, Portugal 
^High Energy Physics, Florida State University, Tallahassee, Florida 32306, USA 
^Instituto de Fisica Corpuscular - IFIC/CSIC, Departamento de Fisica Tedrica 
Universitat de Valencia, 46100 Burjassot, Valencia, Spain 

The simplest unified extension of the Minimal Supersymmetric Standard Model with bi-linear R- 
Parity violation provides a predictive scheme for neutrino masses which can account for the observed 
atmospheric and solar neutrino anomalies in terms of bi-maximal neutrino mixing. The maximality 
of the atmospheric mixing angle arises dynamically, by minimizing the scalar potential, while the 
solar neutrino problem can be accounted for either by large or by small mixing oscillations. One 
neutrino picks up mass by mixing with neutralinos, while the degeneracy and masslessness of the 
other two is lifted only by loop corrections. Despite the smallness of neutrino masses R-parity 
violation is observable at present and future high-energy colliders, providing an unambiguous cross- 
check of the model. 
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The pattern of fcrniion masses and mixings constitutes 
one of the most important issues in modern physics. Here 
we propose a model for the structure of lepton mixing 
which accounts for the atmospheric and solar neutrino 
anomalies. It is based on the simplest one-parameter ex- 
tension of minimal supergravity with bi-linear R-Parity 
violation |l[ as would arise, perhaps, from gravitation. 

The recent announcement of high statistics atmo- 
spheric neutrino data by the SuperKamiokande collab- 
oration has confirmed the deficit of muon neutrinos, 
especially at small zenith angles, opening a new era in 
neutrino physics. Although there may be alternative so- 
lutions of the atmospheric neutrino anomaly |^ it is 
fair to say that the simplest interpretation of the data 
is in terms of to flavour oscillations with maxi- 
mal mixing. This excludes a large mixing among and 
1^, in agreement also with the Chooz reactor data. 
On the other hand the persistent disagreement between 
solar neutrino data and theoretical expectations Q has 
been a long-standing problem in physics. Recent solar 
neutrino data |^ are consistent with both vacuum oscil- 
lations and MSW conversions. In the latter case one can 
have either the large or the small mixing angle solutions, 
with a slight trend towards the latter ||^. The situa- 
tion might become clearer in the near future when rate- 
independent observables such as spectrum, day-night and 
seasonal variations are better measured. In summary one 
sees that while quarks are weakly mixed, there is now 
the intriguing possibility that neutrino mixing is (close 
to) bi-maximal 0. 

Our model breaks lepton number and therefore neces- 
sarily generates non-zero Majorana neutrino masses 
It has strong predictive power and allows for a dynamical 
determination of the atmospheric neutrino angle. More- 
over it leads, under certain circumstances, to bi-maximal 
neutrino mixing. At tree-level only one of the neutrinos 
picks up a mass by mixing with neutralinos [p[, leav- 



ing the other two neutrinos massless While this 

can explain the atmospheric neutrino problem, to recon- 
cile it with the solar neutrino data requires going beyond 
the tree-level approximation. This is the purpose of the 
present paper. For an analysis including only the atmo- 
spheric neutrino problem in the tree-level approximation 
see ref. Q. 

We have performed a full one-loop calculation of the 
neutralino-neutrino mass matrix in the bi-linear 
MSSM. As is shown below, in order to explain the so- 
lar and atmospheric neutrino data it is both necessary 
and sufficient to work at one-loop level. In contrast to 
other papers we have taken special care to achieve 

gauge invariance of the calculation. Moreover we have 
performed the renormalization of the heaviest neutrino, 
thus refining the approximate approaches used, for ex- 
ample, in ref. p^ . For estimates in the approximation 
where loop neutrino masses arise just from tri-linear R- 
parity breaking see ref. p^ ]. 

Bilinear R-parity breaking supersymmetry has been 
extensively discussed in the literature |jl|] . It is motivated 
on the one hand by the fact that it provides an effec- 
tive truncation of models where R-parity breaks spon- 
taneously [|l5| around the weak scale. Moreover, they 
allow for the radiative breaking of R-parity, opening also 
new ways to unify Gauge and Yukawa couplings |lql and 
with a potentially slightly lower prediction for ols [|17| . If 
present at the fundamental level tri-linear breaking of R- 
parity will always imply bi-linear breaking at some level, 
as a result of the renormalization group evolution. In 
contrast, bi-linear breaking may exist in the absence of 
tri-linear, as would be the case if it arises spontaneously. 

Here, we concentrate only on those features of the 
model which are related to neutrino masses. Our model 
consists of the MSSM particle spectrum and superpoten- 
tial except for the addition of the following terms 

W = WMSSM+tiL1Bl. (1) 
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Should supersymmetry not be broken, the above bi-linear 
terms would be superfluous since a suitable redefinition 
of the lepton and Higgs superfields (l8| would convert 
them into trilinear R-parity violating terms. However, 
since supersymmetry must be broken, they give rise to a 
second source for R-parity violation: 



Vsoft = Vsoft.MSSM + BitiL'^H^ 



(2) 



In the presence of soft supersymmetry breaking terms, 
the bi-linear R-parity violating terms can not be rotated 
away, except in the particular case when B = Bi and 
i=l,2,3 which is untypical. Thus we prefer 
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to work in the original basis, containing no trilinear I^p 
vertices. 

The presence of the bi-linear terms in imply that 
the tadpole equations for the sneutrinos are non-trivial, 
i.e. lead to finite VEV for the scalar neutrinos. As a 
consequence the neutrinos and neutralinos, charged lep- 
tons and charginos as well as the Higgses and sleptons of 
the MSSM mix with each other. Detailed mass matrices 
are found in For the neutrino masses the most 

important aspect is, of course, the neutrino-neutralino 
mixing. It generates the following (7 x 7) mass matrix. 



m 



m 



(3) 



where M.^" '^^ usual MSSM neutralino mass matrix 
and the sub-matrix m contains entries from the bi-linear 
parameters. 



(4) 



and (j/g), (t'^) and {vr) are the VEVS of the scalar neu- 
trinos and g' are electroweak gauge couplings. 

It is easy to show that this mass matrix (||) has such a 
structure that only one combination of u^^ i'^, Vt picks up 
a mass, while the remaining two states remain massless. 
This structure is reminiscent of that found in ref. [^0| . If 
the RPV parameters are smaller than the typical size of 
the MSSM parameters, there exists a simple approxima- 
tion formula for the non-zero mass of the neutrino, 



Mig^ + Mag'' 
Adet{M^o) 



|A|^ 



where, 



(5) 



(6) 



and Ml, M2 are supersymmetry breaking electroweak 
gaugino masses. This "alignment" vector plays a promi- 
nent role in all the discussion below since it will fix both 
the overall neutrino mass scale as well as the atmospheric 
neutrino mixing. With two neutrinos being massless one 



of the angles describing the mixing between them can 
be rotated away However, in the presence of loops 
this angle, which will characterize the solar neutrino con- 
versions, will acquire a meaning, together with a (Dirac- 
type) CP phase. 

There are three simple topologies of relevant Feynman 
diagrams contributing to the neutrino-neutralino mass 
matrix . With these the one-loop corrected mass ma- 
trix is calculated as, 



(m) + \ 



(7) 



where and Hy are self-energies. For a complete de- 
scription see . Here, DR signifies the minimal dimen- 
sional reduction subtraction scheme and is the renor- 
malization scale. In order to check for gauge invariance 
in calculating Hy and we have used the general R^ 
gauges. As demonstrated in ref. H^] gauge invariance re- 
quires the inclusion of the tadpole diagrams for the Gold- 
stone bosons associated with the and into the self 
energies. Moreover, in minimizing the scalar potential, 
for consistency reasons it is necessary to also include tad- 
pole diagrams, when solving the tadpole equations, but 
excluding the Goldstone tadpole graphs which have been 
already included into the self energies Q . On the other 
hand, if Goldstone tadpole graphs are kept in the tadpole 
equations rather than in the self energies, gauge depen- 
dent VEVs would be generated. This problem has been 
ignored so far in all previous descriptions | p3| . 

The scalar potential contains terms linear in 
the real part of the neutral scalar fields ctq = 
{a°,a"^,Re{i>i),Re{i}2),Re{i)3)) 



Vlinear = tdC^d + + tiRe{i>i) — ta<7^ 



(8) 



The coefficients of these terms are the tadpoles. Includ- 
ing the one-loop contribution we write 



(9) 



where T^^{Q) = -6t^^ + T^iQ) are the finite one- 
loop tadpoles. The minimization of the scalar potential 
corresponds then to solve ta = 0- This is done by solv- 
ing these equations for the soft masses squared. This 
is easy because those equations are linear on the soft 
masses squared. However the values obtained in this way, 
which we call , are not equal to the values rrif{RGE) 
that we got via the Renormalization Group Equations 
(RGE) starting from universal soft masses at the unifi- 
cation scale. To achieve equality we define a function 



ml{RGE)' 



) 



(10) 
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with the obvious property that 77 > 1. Then we adjust 
the parameters at unification scale to minimize rj. 

We have performed a complete scan of the ^ MSSM 
parameter space, following the procedure outlined above. 
As an example we allow the MSSM parameters to vary 
within the range M2, up to 500 GeV, mo up to 1 TeV, 
and assumed |Ao/too| ^ 3, which helps avoiding charge 
breaking minima. Moreover we assumed tan /3 < 10. The 
latter is needed in order to obtain a nearly maximal atmo- 
spheric angle, since otherwise the sizeable loop involving 
down quarks and squarks would distort this feature due 
to its very strong tan/3 dependence. Note also that the 
bound on tan/3 implies that the lightest CP even Higgs 
boson mass lies below 115 GeV or so. As we will see 
below we can find simultaneous solutions to the atmo- 
spheric and solar neutrino problems only in those parts 
of parameter space where the one-loop contributions to 
the neutrino mass are smaller than the tree-level contri- 
bution. In this case it is possible to give a simple approx- 
imate formula for the composition of the third neutrino 
mass eigenstate. 



Ua.3 — Sin{Atan{ 



-)) 



(11) 
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Figure 1: The atmospheric angle as function of 
|A^/A^|, for \ei\ = e and K = Q.IK- Here e^A < 0.1, 
since larger values lead to larger scatter for very small 
A^/At-]. Maximality of atmospheric mixing is only pos- 
sible for |A^| ~ |A^|. 

Accounting for the atmospheric neutrino anomaly re- 



quires that the Vn 



mixing be large, with little ef- 



fect of Ve in the atmospheric neutrino oscillations. Fit- 
ting for the atmospheric neutrino data then fixes |Ap/A^| 
through this simple equation. These parameters are dy- 
namically determined since they involve Higgs and sneu- 
trino VEVS obtained from the scalar potential. In fact, 
the |Aq,| parameters are simply proportional to the sneu- 
trino vacuum expectation values in the basis where the 



bilinear term in the superpotential is "rotated away" in 
favour of a tri-linear one As an illustration Figure 
1 shows the — Vr angle as a function of |A^/A^| for 
Ae ~ O.lA^ for an otherwise random variation of param- 
eters. Clearly the condition |A^| — \K\ is sufficient to 
ensure near maximal mixing, as long as Ag is somewhat 
suppressed. 

One immediate consequence of the smallness of loop 
with respect to tree contributions is that the absolute 
Ifip scale is then fixed by the atmospheric neutrino mass 
scale. For the above choice of sampling one has |A| ~ 
0.03 — 0.25 GeV^. While this value is surely smaller than 
the weak scale, it may arise naturally in models where 
the sneutrino VEVS are generated radiatively |jl|. 

With the magnitude of parameters fixed by the at- 
mospheric neutrino problem, the question arises, whether 
the loop-induced oscillation parameters, mass splitting 
and angle, are in the right range for either the vacuum or 
the MSW solution to the solar neutrino problem. Since 
the ratio of the loop masses to the tree-level mass de- 
pends on the relative size of the bi-linear Jji^ parameters 
with respect to the alignment vector A this can not in 
general be predicted in the bi-linear model. We have 
found, however, that with our assumption of generation- 
independent bi-linear parameters there should be a 
relative sign between the dynamically determined A 
parameters, i.e. A^ ~ — A.^- Figure 2 shows how, having 
fixed the A^ by the atmospheric neutrino problem, the so- 
lar angle is determined under the above sign assumption 
as a function of Cg/e^. 
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Figure 2: The solar angle as function of ee/^ii, for = 
Cr and Ap = A^, but Ag = O.IA^, applying the condition: 



{A^/Ar) X {ef,/er) < 
only possible for ~ 



0. Maximality of solar mixing is 



As for the solar neutrino scale we show in fig. 3 Amj2 
versus e^/|A|, where = X^i^f- As is seen from the 
figure, for fixed tree-level mass the loop masses depend 
strongly on this quantity. Large values give Am^2 l^e 
MSW range, while low values could give vacuum solu- 
tions to the solar neutrino problem. 
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While a certain amount of "alignment" is needed, for 
masses in the MSW range, the model by itself does not 
prefer one solution over the other. It is also clear that 
bi-maximal neutrino mixing is generated in the bilinear 
$p MSSM - independent of the actual values of SUSY 
parameters, if (i) the I^p bi-linear terms are (nearly) gen- 
eration blind, (ii) tan /9 < 10 implying that the lightest 
CP even Higgs boson mass lies below 115 GeV or so, 
and (iii) ~ — A,-, as long as Ae is somewhat smaller 
than the A^ and A,-. We have checked explicitly that (iii) 
arises dynamically by minimizing the scalar potential of 
the theory. 




eV|A| 

Figure 3: Amfj as function o/e^/|A|. Low values lead 
to neutrino masses in the just-so range, whereas high val- 
ues give Ami2 in the range of the MSW solution. 

Finally, apart from the possible detection of lightest 
CP even Higgs boson (mass below 115 GeV or so), we 
would like to point out that, despite the smallness of neu- 
trino masses R-parity violation is observable at accelera- 
tors through the observation of the decay of the Lightest 
Supersymmetric Particle (LSP), typically a ncutralino. 
For example for a LSP mass of about 50 to 60 GeV the 
decay will occur inside typical LEP, Tevatron and LHC 
detectors (the ncutralino decay length can be of the or- 
der of one meter or so), diluting the missing momentum 
signal to a maximum of 15% of the MSSM expectations. 
The LSP decay will give rise to high multiplicity events, 
providing an unambiguous test of the model. This is 
specially so if branching ratios are measured. In order to 
demonstrate this correlation we have calculated the ratio 
of semi-leptonic branching ratios of the LSP into muons 
and taus. We found the striking result that the correla- 
tion depicted in fig 1 which is required by the atmospheric 
neutrino anomaly is mapped into a well-defined correla- 
tion for the ratio of semi-leptonic LSP branching ratios 
into muons and taus (Figure 4). Note that this correla- 
tion holds for the semi-leptonic decays, despite the fact 
that there are many scalar boson exchanges contributing 
to the LSP decay! The case where the LSP is heavier 
than the W and has 2-body decays mainly to W and Z 



was considered in ref. who found a similar correlation 
in the approximation where the LSP decay into lightest 
neutral supersymmetric Higgs boson was neglected. The 
result we show in Figure 4 is general, independently of the 
ncutralino mass, and provides a powerful way to probe 
the solution of the atmospheric neutrino anomaly and 
opening the potential to measure the related neutrino 
angles at high energy accelerators! 
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Figure 4: Ratio of branching ratios for semilep- 
tonic LSP decays into muons and taus: BR{x —> 
liq'q)lBR{x — > Tq'q) as function o/ A^/A.^- 
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